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' , Abstract. We compare recent lattice results on fluctuations and correlations of strangeness, 

ijT^. baryon number and electric charge obtained with p4 improved staggered action with the 

(-H ' prediction of hadron resonance gas model. We show that hadron resonance gas can describe 

these fluctuations reasonably well if the hadron properties are as calculated on the lattice. 

> 

^■r^ , 1. Introduction 

^D ' In recent years a lot of progress has been made in understanding the QCD transition at non-zero 

l/-^ ■ temperature through lattice calculations (see Refs. [HE] for reviews). In particular, fluctuations 

^D \ of conserved charges have been studied on the lattice and provided insight on how the relevant 

degrees of freedom change from hadronic to partonic (Sj HJ (Sj (U [71 [U |9l [10] . At low temperatures 
it is expected that hadrons are good degrees of freedom and thermodynamic quantities can be 
well described by hadron resonance gas. Hadron resonance gas (HRG) turned out to be very 
^ ■ successful in describing particle abundances produced in heavy ion collisions [H]. It was also 

used to estimate QCD transport coefficients [12] as well as chemical equilibration rates [13] 
close to the transition temperature. Thermodynamic quantities calculated in lattice QCD with 
rather large quark mass agree well with the HRG model if the masses of the hadrons in the 
model are tuned appropriately to match the large quark mass used in lattice calculations |14j . 
Furthermore, the ratios of certain fluctuations calculated on the lattice are also in reasonably 
good agreement with HRG model predictions at low temperatures |1 ^ [T6 | [TO] . 

Fluctuations of conserved charges are the most suitable quantities to test the validity of HRG 
model. They can be defined as derivatives of the pressure with respect to chemical potentials 
and as such appear in the calculations of thermodynamic quantities at finite baryon density 
via Taylor expansion. The recent lattice calculations performed with p4 and asqtad actions at 
or very close to the physical quark mass [6l [3 [TO] gave fluctuations that were quite different 
from the results obtained in HRG model. It was pointed out that lattice discretization effects 
on the hadron spectrum are responsible for this discrepancy [1^. When taking into account 
the lattice spacing dependence of the hadron masses in the hadron resonance gas calculations a 
good agreement between lattice and HRG calculations has been found. Strangeness fluctuations 
calculated with and stout [8l [9] and HISQ [18J improved staggered fermion actions agree better 
with the HRG result (see also the discussion in Ref. [IS]). This is due to largely reduced 
discretization effects in the hadron spectrum for these actions. In Ref. [IT] we studied baryon 
number and strangeness fluctuations in HRG, where the hadron masses have been modified to 
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include discretization effects present in lattice calculations. The comparison was done using 
lattice data obtained with asqtad action on A'^^ = 6 and 8 lattices. On the other hand, the 
most detailed lattice study of baryon number, strangeness and electric charge fluctuations and 
correlations up to sixth order was performed with p4 action using lattices with temporal extent 
Nr = 4 and N-j- = 6 [TO]. In this contribution we are going to extend our calculations in 
modified the HRG model to study different fluctuations and compare them with lattice results 
obtained with p4 action on A',- = 4 and 6 lattices. As mentioned above this is also important 
for constructing realistic equation of state at finite baryon density for hydrodynamic models, 
similarly as this was done for zero baryon density in Ref. [17j . 

2. Fluctuations of strangeness, baryon number and electric charge 

Derivatives of the pressure with respect to chemical potentials of conserved charges, e.g. baryon 
number {B), electric charge Q and strangeness S at zero chemical potentials can be well 
calculated in lattice QCD 

..X _rT.n d''p{T,HB,t^Q,Hs) ^ y _ T, n Q (^\ 

These are related to quadratic and higher order fluctuations of conserved charges X2 — 
(X2)/(yr3), xf = ((A^x) - 3(A^i)^)/(^T3) etcEl Mixed derivatives of the pressure give 
correlations of conserved charges 

We have calculated these quantities in HRG with hadron spectrum modified to take into account 
the lattice artifacts. Unfortunately, there is no detailed calculation of the hadron spectrum with 
p4 action. Therefore the lattice spacing and quark mass dependence of the hadron masses was 
evaluated using the formulas derived in Ref. [T7] based on asqtad calculations. We expect that 
cutoff dependence of the hadron masses is similar for p4 and asqtad actions, except in the pseudo- 
scalar meson sector. It is known that the quadratic splitting of non-Goldstone pseudo-scalar 
mesons is about two times larger for the p4 action than for the asqtad action ^20j. Therefore 
when evaluating the contribution of kaons and pions to different quantities we simply doubled 
the quadratic splittings in the pseudo-scalar meson sector used in Ref. [T7]. To get agreement 
between lattice results and HRG it turned out to be necessary to modify the masses of excited 
baryons states. Since the cutoff dependence of the excited baryon masses is not known, it was 
assumed that masses of all excited states up to certain threshold m^^ are modified the same 
way as the ground state baryon masses yL7j. Values m^^ = (1.7 — 2.5)GeV have been considered 
in the previous analysis. 

In this paper we consider all resonances with mass up to 2.5GeV and use m^^^ = 1.9GeV 
in all calculations. The quadratic fluctuations of strangeness and electric charge calculated in 
HRG and compared to the p4 lattice results are shown in Fig. [TJ The lattice results are well 
below the HRG curve obtained with physical hadron masses. This is a general feature for all 
fluctuations. However, taking into account the discretization effects in the hadron spectrum we 
flnd a reasonable agreement between HRG and lattice results. In the strange sector discretization 
effects are slightly overestimated in our approach. We also studied quadratic baryon number 
fluctuations and baryon number - strangeness correlations and the corresponding results are 
shown in Fig. [2l The baryon number fluctuations are well described by modified HRG, 
although the in the low temperature region cutoff effects are under-predicted for Nj- = 6. For 

^ Here we consider the case of zero chemical potential, so (A'^x) = 0. 
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Figure 1. The quadratic fluctuations of strangeness (left) and electric charge (right) calculated 
on the lattice with p4 action pi)j and compared with the HRG calculations in the continuum 
(solid black line) and on the lattice (blue and red lines). 
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Figure 2. The quadratic fluctuations of baryon number (left) and baryon number - strangeness 
correlations (right) calculated on the lattice with p4 action [10] and compared with the HRG 
calculations in the continuum (solid black line) and on the lattice (blue and red lines). 

baryon number-strangeness fluctuations the agreement between lattice and modified HRG is 
not very good. This may imply that the cutoff effects in strange baryon sector are smaller than 
in non-strange baryon sector and are overestimated in our approach. Interestingly, the new 
calculations with HISQ action also suggest smaller cutoff effects in baryon number-strangeness 
fluctuations compared to baryon number fluctuations p^. Finally we have considered fourth 
order fluctuations of the electric charge and baryon number. The numerical results are shown 
in Fig. [3l The HRG calculations with the modified hadron masses agree reasonably well with 
the p4 lattice data. Note, however, that here the deviations from the resonance gas show up at 
smaller temperatures. 



3. Conclusions 

We have studied quadratic and fourth order fiuctuations of baryon number, electric charge and 
strangeness fiuctuations in HRG model, where hadron masses have been adjusted to take into 
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Figure 3. The fourth order fluctuations of the electric charge (left) and baryon number (right) 
calculated on the lattice with p4 action [TO] and compared with the HRG calculations in the 
continuum (solid black line) and on the lattice (blue and red lines). 

account the discretization errors in the hadron spectrum present in the lattice calculations. We 
have found reasonably good agreement with the lattice calculations of the fluctuations performed 
using p4 action on Nr = 4 and 6 lattices. Our calculations explain why all the lattice results 
fall below the physical HRG result. Our approach does not give a good description of baryon 
number - strangeness fluctuations. To resolved this issue a more refined treatment of the cutoff 
effects in the baryon sector is needed. 
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